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The modulated beam of the free electron laser when 
passed through a thin target should produce strong 
transition radiation. The measurement of the transition 
radiation is a direct measurement of the electron beam 
modulation in the free electron laser. A transition 
radiation experiment using the Stanford MKIII Infrared Free 
Electron Laser (IRFEL) has been proposed. The analysis has 
centered on TRANSPORT, a computer program used for designing 
charged particle beam systems. The MKIII IRFEL wiggler exit 
bending magnet system was modeled using TRANSPORT. Analysis 
reveals that the transverse emittance and momentum spread 
characteristics will cause the modulated beam to demodulate 
along the path of the central trajectory. A detector 
location 10 millimeters downstream of the first bending 
magnet is found suitable for the measurement of the electron 
beam modulation. For this case the thin foil must be 
rotated approximately 11.25 degrees about the y-axis to 
yield a minimum effective picobunch extent as seen by the 
thin foil. Alternate FEL systems are suggested for 
investigation. 
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I. 



INTRODUCTION 



A. OVERVIEW 

The electron beam from a free electron laser (FEL) 
offers the possibility of observing and using transition 
radiation (TR) at optical frequencies. After interacting 
with the wiggler and radiation fields, the electron beam 
becomes modulated in the axial direction at the lasing 
frequency [Ref. 1]. The axially modulated bunches 
(picobunches) could be passed through a thin foil producing 
the surface effect TR as illustrated in Figure 1-1. 



back\ 

specu 



beam 

picob 




forward TR 



viewport 



Figure 1-1 Transition Radiation (TR) 
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The axially modulated beam picobunches impact a thin foil of 
aluminum or mylar yielding forward and backward/specular TR 
[Ref. 2]. Both forms of radiation should be enhanced by 
coherent emission from all charges in the picobunch. The 
backward/specular TR is redirected to a viewing port/optical 
detector by a mirror. The resulting radiation could be used 
to diagnose the electron beam modulation [Ref. 1] . 

The TR power radiated per picobunch is found to be 
proportional to the charge of the picobunch squared. For a 
picobunch of approximately 10 6 electrons, the radiation 
intensity can easily be detected and measured with an 
optical detector. The axially modulated electron beam of 
the IRFEL, upon impact with the thin foil or multifoil, will 
generate coherent optical transition radiation. [Ref. 2] 

B. MODULATION OF THE ELECTRON BEAM 

The longitudinal bunching of the electron beam by the 
undulator of the free electron laser (FEL) is a result of 
the force acting on the individual particles by the electric 
field of the superimposed laser radiation. This effect on a 
single electron is examined in Figure 1-2 [Ref. 3]. This 
figure illustrates the fields and trajectories within the 
undulator of a FEL: an alternating magnetic field, a co- 

propagating laser electric field and an electron along a 
trajectory. The electric field of the laser is zero at 
position A. The electron experiences no force other than 
that of the undulator' s magnetic field, which will bend it 
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Figure 1-2 Longitudinal Bunching of the Electron Beam 



in the horizontal plane. Now consider the electron a time 
later at position B. The electron has traveled one-quarter 
of an undulator period and now experiences a maximum in the 
laser electric field. The force on the electron at point B 
is opposite the direction of the laser electric field at 
that point and is labeled F e . The electron will therefore 
experience a decelerating force. Whether an electron is 



accelerated 


or 


decelerated by 


the 


laser electric 


field 


depends 


on 


the 


phase of the 


laser 


electric field 


with 


respect 


to 


the 


electron oscillation. 


Some electrons are 



accelerated, while those half a laser wavelength ahead or 
behind are decelerated. For example, an electron at C would 
be accelerated. As the accelerated electrons catch up to 
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those which have been decelerated, they form into bunches at 
the laser wavelength. In other words, the electron beam 
becomes axially modulated at the laser wavelength of the FEL 
radiation. [Ref. 3] 

Figure 1-3 [Ref. 4] depicts approximately 10 6 electrons 
being bunched within the optical wavelength of the FEL 
radiation while within the undulator. 




ELECTRONS 



UNDULATOR 




| 4 X ► 




BUNCHED ELECTRONS 



WAVELENGTH OF LASER RADIATION 



Figure 1-3 Bunching of Electrons within the Wavelength 
of Laser Radiation 

As the electrons oscillate they spontaneously emit "magnetic 
bremsstrahlung" peaked in the forward direction thus 
generating the laser electric field. When the laser 
frequency is nearly resonant with the electron oscillations, 
significant amounts of energy can be exchanged between the 
electrons and the laser field causing the bunching effect of 
the electron beam. The condition for this resonance is 
given by Equation 1-1 (MKS) [Ref. 5]: 
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( 1 - 1 ) 



Y 2 — A W /2X(1 + A^^/4 Tr^itiQ^Q^ ) , 



where: 

y = energy of resonant electron divided by rest 
mass energy, 

> w = wiggler (undulator) magnetic period, 

X = laser wavelength (operating wavelength) , 
e = electron charge, 

Bj-ms = rms magnetic induction of the wiggler, 
c = velocity of light in vacuo. 

The operating (lasing) frequency of a FEL is determined 
by Equation 1-2 (MKS) , 

X = X w /2y 2 (1 + K 2 ) , (1-2) 



where: 



K = the dimensionless measure of the wiggler 
strength , 

K = eB rinS A w /2'rrm 0 c . (1 — 3) 

The operating parameters of the Stanford MKIII IRFEL are 
shown in Table 1-1 [Ref. 6]. 

The wiggler period of the IRFEL is of the "hybrid" type 
employing both SmCo permanent magnets and vanadium permendur 
steel poles. The magnets are oriented to provide a nearly 
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TABLE 1-1 



MKIII IRFEL PARAMETERS 



WIGGLER 








magnet period 


A w 


2 . 3 


cm 


number of periods 


N 


47 . 5 




minimum gap 


9min 


6.5 


mm 


wiggler strength 


K 


1.14 




peak magnetic field 


Bo 


0.702 


T 


RESONATOR 








length 


L 


178 . 57 


cm 


Rayleigh range 


Z R 


73 


cm 


E-BEAM 








energy 


Y 


85 




energy 


E 


43 


MeV 


energy spread 


sy/y 


0.5 


% 


peak current 


I 


20 


Amps 


emittance 


e N 


7 


Tr-mm mrad 


macropulse length 


T 


3 


lisec 



sinusoidal vertical field By, using a Halbach undulator 
design. [Ref. 6] The lasing frequency can be modified by 
slightly changing the energy of the beam or by changing the 
K value of the wiggler. In practice the K parameter is 
modified by adjusting the gap between the permanent magnets 
in the wiggler cavity. The resonator length must then be 
modified to account for the new wavelength. Assuming the 
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lasing wavelength of 3.0 microns, the IRFEL will produce a 
modulated beam as shown in Figure 1-4 . 




3u 





Figure 1-4 Disks of Charge (picobunches) 

The rectangles indicate disks of charge of thickness 1 
micron separated by 3 microns. This periodic bunched pulse 
train will generate transition radiation of wavelength 3 
microns and also odd harmonics when passing through a thin 
foil such as mylar or aluminum. The current vs time profile 
is shown in Figure 1-5. 
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C. WHAT IS TRANSITION RADIATION 

Transition radiation (TR) is the electromagnetic radiation 
which is emitted when a uniformly moving charged particle 
traverses from one medium into another with a different 
dielectric constant. More generally the effect will occur 
in the presence of a heterogeneity in a medium. [Ref. 7] 

TR may be considered to originate from the change in the 
fields as an electron crosses a boundary. The magnetic 
field remains essentially constant but the electric field 
changes because of the change in the dielectric constant. A 
third field, radiation, is required to meet the boundary 
conditions . 

TR can also be explained in terms of a charge and its 
image. While an electron moves uniformly in vacuo it does 
not radiate; the field in vacuo is equal to the field of the 
electron and its image moving toward it. When the electron 
crosses the boundary surface both the electron and its image 
cease to exist from the point of view of the field. The 
crossing of the boundary must give rise to exactly the same 
radiation as that which is due to a sudden stop at the same 
point of the electron and its image. [Ref. 7] 

The free electron laser will longitudinally modulate the 
electron beam in the undulator/wiggler yielding disks of 
charge of small thickness (lu) . Due to the low efficiency 
of OTR, 10^ electrons per picobunch will yield approximately 
10 4 photons per picobunch upon foil impact. 

For a relativistic beam of electrons ( $ ->■ 1) , and a TR 
foil inclined at 45 degrees with respect to the beam axis 



8 



(Figure 1-1) , the backward/specular reflection intensity for 
a vacuo/dielectric interface is given by the relation [Ref. 
2 ] : 



d 2 I/dwd Q = F( ip,w) (e 2 /tr 2 c) [6 2 /(y -2 + 6 2 ) 2 ] , (1-4) 



where : 



r [ e-sin 2 (ip+0 )] - e cos (ip+6 )_ ■, 2 

f ( ([)/ w) = I 2 T72 T J 

[e-sin (i^+6)J 7 + e cos(<J>+6) 



F = Fresnel coefficient for reflection of light 
with the E-vector in the plane of incidence, 

\p = angle of specular reflection (45 degrees) , 

e = dielectric constant of foil material. 



The intensity of this radiation should be sufficient for 
the detection and observance of TR from a single foil/multi- 
electron interface. Measurement of the radiation will be 
accomplished via a detector of the phototube, TV camera, or 
reticon type. 

The characteristics of TR that will apply to the 
proposed experiment are: 

(1) The angular distribution of TR is such that for 
relativistic particles the emission angle of the 
photons with respect to the electron beam is 
appproximately equal to 1/y . 

(2) The intensity of TR is a strong function of the 
energy of the particle producing it. 
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( 3 ) 



The radiation is polarized with the electric field in 
the plane containing the particle velocity vector and 
the wave propagation vector. 

(4) In the case of a periodic medium (i.e., multifoil) an 
rir phase difference is the condition for constructive 
interference or resonance between radiation produced 
at successive interfaces; r is a positive integer. 
For m boundaries the intensity varies as m 2 . 

(5) If the medium is a transparent dielectric then 
e'=n 2 (w), where n is the refractive index and e' is 
the dielectric constant. 

(6) For a relativistic electron, TR will arise at the 
surface of the interface. [Ref. 7] 

(7) The power radiated from a transition foil is 
proportional to the charge squared [Ref. 2]. 
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II. PROPOSED EXPERIMENT USING MKIII IRFEL 



A. EXPERIMENT PROPOSAL 

The purpose of this experiment will be to measure the 
electron beam modulation of the FEL using a TR foil. 

The MKIII IRFEL, located at the High Energy Physics Lab 
at Stanford University may lend itself to an experiment to 
measure the axial modulation of the electron beam of the 
FEL. The "bunched" beam downstream of the wiggler magnet is 
normally bent away from the downstream cavity mirror and 
directed to the beam dump. If the beam were to retain its 
axial modulation, upon the passage through a thin foil it 
would generate coherent transition radiation. The thin foil 
will hereafter be referred to as the detector. The detector 
must be conveniently located downstream of the wiggler. 

1 . System Description 

The MKIII IRFEL is shown in Figures 2-1, 2-2, 2-3, 

and 2-4. Figure 2-1 is a photo of the accelerator section 
of the IRFEL. Figure 2-2 is a representation of the IRFEL 
downstream of the accelerator. The electron beam produced 
by the accelerator is directed around the upstream optical 
cavity mirror and through the wiggler magnet. The beam is 
axially modulated within the wiggler and is then directed 
around the downstream movable mirror to the beam dump. An 
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Figure 2-1 MKIII IRFEL Accelerator 



12 





nr 

vv. 

.*•/ 






• s*. 









■V,l 






















rW 

w: 






yr»! 




















,-V 

V*, 

■>V; 






$ 

£ 





13 



Figure 2-2 MKIII IRFEL 





external photo of the MKIII IRFEL wiggler magnet is shown in 
Figure 2-3. 




Figure 2-3 MKIII IRFEL Wiggler 



Figure 2-4 is a photo of the experimental beamline apparatus 
downstream of the wiggler cavity. The electron beam will 
proceed from left to right. The figure shows the waveguide 
type piping, flanges, bending magnet #2, bending magnet #3, 
and a four way cross. Bending magnet #1 is upstream to the 
left of the flange and bellows (not shown in the photo) . 
The detector could easily be placed on the waveguide piping 
or at the four way cross location. 
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Figure 2-4 MKIII IRFEL Wiggler Exit 

2 . Possible Detector Locations 

Figure 2-5 is a top view representation depicting 
possible locations for the detector. The electron bunches 
exiting from the wiggler are deflected by the first bending 
magnet in order to prevent them from striking the resonator 
mirror. 

3 . Hardware Design at Detector Locations 

The hardware design needed for a detector between 
locations B and C is shown in Figure 2-6. Placement of the 
detector between points B and C would require modifying the 
waveguide piping, bellows, and flange. Location E is the 
alternative location of the detector and would require the 
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Figure 2-5 Detector Locations 




Figure 2-6 Detector Placement on Waveguide 
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modification of Figure 2-7A. The linear motion feedthrough 
needs to be connected to a sliding sleeve device in order to 
position the TR foil and mirror. A suitable sliding sleeve 
device is shown in Figure 2-7B. In general, it would be 
easier to modify the IRFEL by placing the detector at 
location E. 

B. ANALYSIS OF ELECTRON BEAM TRANSPORT EQUATIONS 

Note: Readers not familiar with the properties of beams and 
the TRANSPORT development of system equations should read 
Appendices A and B before proceeding with this section. 

The development of the electron beam transport equations 
for locations A through E along the central trajectory is 
found in Appendix B. It is desired to find the best 
location in the beamline to place the detector in order to 
generate coherent transition radiation. Of particular 
concern is the longitudinal extent of the picobunches 
(Figure 1-3) as they traverse the beamline. In the ideal 
case the picobunches would retain their axial modulation 
traversing the beamline. This is the case of a non- 
divergent, monoenergetic beam. In the real world, electron 
beams are non-monoenergetic and are divergent. Beams of 
this nature are specified in terms of their phase space. In 
order to observe coherent transition radiation, the 
picobunches must not "debunch" significantly upon traversing 
the beamline. "Significant" debunching for this experiment 
is on the order of a few (>3) microns. If debunching 
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Linear motion 
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Beam in 



Viewport 
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Figure 2-7A Detector Placement at Location E 
(All Dimensions in Inches) 
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Figure 2-7B Sliding Sleeve Device 

(All Dimensions in Inches) 
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exceeds 3 microns the coherent TR effect will not be 
observed. A simple representation of the beam traversing 
the system is desired. Figure 2-8 is a representation of 
the paths a single electron located at "a" may take in 
traversing the system. System parameters (obtained from 
blueprint) and definitions applicable to the representation 
are shown below the figure. 




Figure 2-8 Trajectory of a Single Electron 



x 0 = the initial horizontal displacement of the 
electron with respect to the central trajectory. 

6 o = the divergence of the electron with respect to 
the assumed central trajectory. 

<5 0 = the fractional momentum deviation of the 
electron from the assumed central trajectory. 
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The central trajectory parameters are: 



= drift length 1 
L 2 = drift length 2 



327.03 mm 



95.30 mm 



L = effective pathlength through 
bending magnet 



38.35 mm 



a = bending angle 



11.25 0 



Po = bending radius of magnet 
L 0 = initial drift length 
E = energy of electron 



195.30 mm 



0.0 mm 



40 MeV 



(Note: L 0 = 0; the electron starts at the magnet pole face) 

1 . Buskirk System Equations 

The system equations for this model were derived 
by Dr Fred Buskirk. These equations assume the 
initial drift length, L 0 and the effective 
pathlength through the bending magnet, L, equal 
zero. The derived equations are: 

dx = (Li/cos 2 a + L 2 )0 o + (L^sin a/cos 2 a) S Q + x Q , (2-1) 
ds = - (L 1 sina/cos 2 a)0 Q - (L 1 tan 2 a)6 0 , (2-2) 



dx = the change in the horizontal distance of the 
electron with respect to the central trajectory, 
and 

ds = the difference in the pathlength between the 
arbitrary ray and the central trajectory. 



where : 
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The electron is to traverse the system from point "a" to the 
xy plane of location E. Assuming the fringing fields are 
negligable the electron will bend in the xz plane only. 

Trajectory 1 (Tl) is for the case of the 
monoenergetic ( 6 0 = 0) non-divergent (0 O = 0) electron 
initially located at x 0 . Its trajectory coincides with the 
central trajectory as expected. Trajectory 2 (T2) is the 
case of a non-monoenergetic non-divergent electron initially 
located at x Q . Since this electron possesses greater 
momentum than the central momentum it will not be bent as 
much by the bending magnets . The pathlength between the 
bending magnets will be shorter than . If two electrons 
having equal velocities were located at "a" start downstream 
on trajectories Tl and T2 , the electron on T2 would have a 
shorter pathlength through the system. At the time when the 
electron on Tl reached the xz plane of location E, the 
electron on T2 would have already passed throught the system 
by a distance equal to the pathlength difference. T2 is 
also slightly greater than Tl in the x direction. 
Trajectory 3 (T3) is the case of a non-monoenergetic, 
divergent electron initially located at x Q . As this 
electron traverses the system on T3 , the pathlength between 
the bending magnets will be different from the central 
trajectory and the electron would end up at a different 
location downstream. The three electron trajectories 
illustrate the change in pathlength due to the effects of 
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divergence and percent momentum change. This is analogous 
to the change in the longitudinal extent of a picobunch as 
it traverses the system. 

2 . TRANSPORT Equations 

The TRANSPORT equations for location E are derived 
in Appendix B. These equations are: 

x E (t) = (2p 0 tana+ L^cos 2 a + L 2 ) 0 O 

+ ( -2P 0 + 2 p^/cos a + L 1 sinc/cos 2 a)6 Q , (M-13) 

l E (t) = ( — 2p 0 + 2 f^/cosa +L!sin a/cos 2 a) 0 Q 

+ (- 2 p 0 a + 2 p o tan 0 t + L^tan 2 a)6 0 + lo* (M— 15) 

It can be shown that the TRANSPORT equations reduce to the 
Buskirk equations letting the effective pathlength through 
the bending magnet, L, and the initial one-half longitudinal 
extent of the picobunch, 1 Q , equal zero. The variables x 
and 1 correspond to dx and ds respectively in Equations 2-1 
and 2-2. TRANSPORT analysis requires the phase space 

representation of the input beam. For a normalized 
horizontal beam emittance of 8 mm-mrad (revised MKIII IRFEL 
parameter) with a beam energy of 40 MeV, assume the 2-D 
phase space of Figure 2-9. The points a-d corresponding to 
the extreme values of the phase space ellipse are to be 
plotted. This will yield a physical representation of what 
the detector will see at the location of interest. 
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Substituting the system parameters with 1 Q = .0005, 
Equations M-13 and M-15 become: 

X ( t ) = .512960 O + . 739766 Q + x Q , (2-3) 

l(t) = .O73980 o + . 13940<5 o + .0005 , (2-4) 

where x Q/ x(t), 1 0/ l(t) are in mm; 0 O in mrad, 6 Q in 
% A P/P. The plot of x(t) vs l(t) for the points in Figure 
2-9 is shown in Figure 2-10. Tabulated values are shown in 
Table 2-2. The center ellipse corresponds to a mono- 
energetic picobunch with divergence in the xz plane. The 
divergence alone causes a significant increase in the 
picobunch extent. A picobunch of longitudinal extent 1 
micron grows to a longitudinal extent of approximately 14 
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Fiaure 2-10 Picobunch Extent at Location E 



microns at location E. The right side ellipse represents 
the same picobunch having one-tenth of one-percent greater 
momentum and the left side ellipse represents the picobunch 
with one-tenth of one-percent less momentum. A non- 
monoenergetic picobunch with 6 0 = +/- .1% (revised MKIII 

IRFEL parameter) would have a total longitudinal extent of 
= 4 6 microns at location E. A coherent TR effect would not 
be expected under these conditions. Even in the mono- 
energetic case, the beam divergence would have to be 
greatly reduced to yield a longitudinal extent of <= 3 
microns. The TRANSPORT system equations for detector 
placement midway between locations B and C are: 

Xc(t) = (cosa)x 0 + (p Q sina + Li/cosa)0 o 

+ [p Q (l-cosa) + Litana]<5 Q , (M-5) 

left) = (-sim)x 0 - p 0 ( 1-cos a) 0 O + 1 Q 

- p 0 (a-sina)S 0 . (M-7) 

Substituting the known parameters with = 163.52 mm these 
equations become: 

X ( t ) = . 98 079x o +. 204820 0 + .36278 ^ , (2-5) 

l(t) = - . 19509x o - .00375 0 O -.00246 6 o + .0005 . (2-6) 
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A plot of l(t) vs x(t) for the phase space of Figure 2-9 is 
shown in Figure 2-11. Tabulated values are shown in Table 

2-3 . 

Inspection of Figure 2-11 reveals a shift in the 
orientation of the ellipses with respect to Figure 2-10. An 
11.25 degree rotation of the thin foil about the y-axis will 
minimize the total longitudinal extent of the original 
picobunch as seen by the thin foil. The effective 

longitudinal extent of the original picobunch is 

approximately 23 microns for constant x(t) . Moving the 
detector further upstream to a location 10 mm downstream of 
the first bending magnet the TRANSPORT system equations are: 

X ( t ) = . 98079x o + . 04830 0 O + .05742<5 Q , (2-7) 

l(t) = - . 19509x o - . 00375 0 Q - .00246<$ o +.0005 . (2-8) 

A plot of l(t) vs x(t) for this detector location is shown 
in Figure 2-12. Tabulated values are shown in Table 2-4. 
Note that Equation 2-6 and Equation 2-8 are identical, thus 
the longitudinal extent of the picobunch is constant between 
B and C. However, the effective picobunch extent will 
increase between B and C. Figure 2-12 reveals that a 
momentum spread of one-tenth of one-percent will not cause 
significant increase in the effective longitudinal extent of 
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Figure 2-11 Picobunch Extent Halfway Between Bending Magnets 
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Figure 2-12 Picobunch Extent at 10 nun 



the picobunch. The non-monoenergetic ellipses appear to 
closely coincide with the central momentum ellipse. 

It is apparent from equations M-7 and M-15 that the 
change in the longitudinal extent of the original picobunch 
is due to the initial horizontal displacement of the 
electron, x Q , the divergence of the electron, 0 O , and the 
momentum increment, <5 Q ; all with respect to the assumed 
central trajectory. A scaled plot is used to illustrate the 
relative effect of each factor contributing to the change in 
the horizontal and longitudinal extent of the picobunch as 
it traverses the system. Substituting the known parameters, 
the TRANSPORT equations to be plotted reduce to the 
equations of Table 2-1. 



TABLE 2-1 

EQUATIONS FOR TRAJECTORY PARAMETERS 
AT VARIOUS LOCATIONS 

Detector 10 mm from bending magnet #1: 

X ( t ) = . 98079x o + .048309 o + .05742 
l(t) = - . 19509x o - .003756 o + 1 0 “ . 00246 6 0 
Detector halfway between bending magnets #1 and #2: 
X ( t ) = . 98 07 9x 0 + .2O4820 o + .362786 0 
l(t) = - . 19509x o -.OO3750 o + 1 Q -.00246S o 
Detector at location E: 

X ( t ) = X Q + .5129600 + . 739766 q 
l(t) = .0739800 + 1 Q + .13940S o 
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Figure 2-13 is an unsealed representation of the 
effect of the initial horizontal displacement, x 0 , on the 
longitudinal extent of the beam. The beam parameters are 
specified by the given 2-D phase space. Figure 2-13 
illustrates that the initial horizontal displacement with 
respect to the central trajectory will have an effect on the 
longitudinal extent of the beam between the magnets but will 
have no effect beyond the second magnet at location E. 

Figure 2-14 is a scaled representation of the effect 
of the divergence, 0 O on x(t) and l(t). The effect of the 
divergence term on l(t) is small over the regions between 
the magnets but increases beyond the second bending magnet. 
Note that the change in the longitudinal extent is measured 
with respect to the beam direction. 

Figure 2-15 illustrates the effect of the momentum 
increment, S 0 , on x(t) and l(t). 

Comparing Figures 2-14 and 2-15, it appears that the 
divergence and momentum increment have roughly the same 
effect in debunching the beam picobunches as they traverse 
the system. The previous three figures and Figures 2-10, 2- 
11, and 2-12 lead to the physical representation of the 
beam picobunches shown in Figure 2-16. 
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Figure 2-13 Effect of 
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Figure 2-14 Effect of 0 on x(t) and l(t) 
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Figure 2-15 Effect of 
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Figure 2-16 Physical Representation of Beam Picobunch 



TABLE 2-2 



BEAM PARAEMTERS FOR DETECTOR AT LOCATION E 



Point 


(x, 6) 


% Ap/p 


1ft) mm 


x ( t ) mm 


a 


(1,0) 


0 


.001 


1.000 






. 1 


. 015 


1.074 






-.1 


-.013 


.926 


b 


(0, .1) 


0 


. 008 


. 051 






. 1 


. 022 


. 125 






-.1 


-.006 


-.023 


c 


(-1,0) 


0 


. 001 


- 1.000 






. 1 


. 015 


-.926 






-.1 


-.013 


-1.074 


d 


(0,-.l) 


0 


-.006 


-.051 






. 1 


. 008 


. 023 






-.1 


-.020 


-.125 





TABLE 2-3 






BEAM 


PARAMETERS 


FOR DETECTOR HALFWAY 






BETWEEN BENDING 1 


MAGNETS 




Point 


(x,„ 9) 


% AP/P 


1 ft) mm 


x f t ) mm 


a 


(1,0) 


0 


-.194 


. 981 






. 1 


-.194 


1.017 






-.1 


-.194 


. 945 


b 


(0, .1) 


0 


. 001 


. 020 






. 1 


0 


. 057 






-.1 


. 001 


-.016 


c 


(-1,0) 


0 


. 196 


-.981 






. 1 


. 196 


-.945 






-.1 


. 196 


-1.017 


d 


O 

1 

• 

H 


0 


. 001 


-.020 






. 1 


. 001 


. 016 






-.1 


. 002 


-.057 
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TABLE 2-4 



BEAM PARAMETERS FOR DETECTOR 10 mm 
DOWNSTREAM OF FIRST BENDING MAGNET 



Point 

a 



b 



c 



d 



(x , 9 ) 

( 1 , 0 ) 

( 0 , . 1 ) 

(- 1 , 0 ) 

( 0 / ” • 1 ) 



% AP/P 
0 

. 1 
1 

0 

.1 

1 

0 

. 1 
-. 1 

0 

. 1 
-. 1 ' 



l(t) mm 

-.194 

-.194 

-.194 

. 001 
0 

. 001 

. 196 
. 196 
. 196 

. 001 
. 001 
. 002 



x ( t ) mm 

.981 

.987 

.975 

. 005 
. Oil 
-.001 

-.981 

-.975 

-.987 

-.005 
. 001 
-.011 
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III. RESULTS AND CONCLUSIONS 



The axially modulated beam of the MKIII IRFEL was 
represented as disks of charge (picobunches) as shown in 
Figure 1-4. If these picobunches were to retain their axial 
modulation and pass through a thin foil, coherent optical 
transition radiation would be generated. 

The analysis of the TRANSPORT system equations for the 
detector locations of Figure 2-5 reveal that significant 
debunching (of more than 3 microns in the axial direction) 
will occur at these detector locations due to the effects of 
a non-monoenergetic, divergent beam. 

For detector placement at the four-way cross location 
(location E of Figure 2-5) the increase in the picobunch 
extent is approximately 46 microns (Figure 2-10) . For the 
detector location halfway between bending magnets #1 and #2 
(Figure 2-5) the effective picobunch extent is reduced to 23 
microns (Figure 2-11) with the detector rotated 11.25 
degrees about the y-axis. Placing the detector further 
upstream has the overall effect of reducing the effective 
picobunch extent. The axial modulation of the electron beam 
of the MKIII IRFEL can be measured at a detector location 10 
millimeters downstream of the first bending magnet with the 
detector rotated 11.25 degrees about the y-axis. At this 
location there will be a 2 micron increase in the effective 
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picobunch extent (Figure 2-12) . In theory, detector 
placement at this location will enable the desired OTR 
effect to be observed. This location, however, is not 
easily accessible due to structural hardware upstream of the 
bellows of Figure 2-4. In addition, the bellows location 
and waveguide material would have to be modified in order to 
mount a detector of the type illustrated in Figure 2-6. To 

-v*’ 

situate a detector at this location it will be necessary to 
make hardware modifications to accomodate the experimental 
arrangement. 

The normalized horizontal and vertical emittance 
parameters for this experimental proposal were 8 mm-mrad and 
4mm-mrad respectively. The momentum spread was one-tenth of 
one percent. A reduction in these parameters will yield a 
reduced picobunch extent at any detector location. 

In addition to the Stanford MKIII IRFEL, there are 
several IRFELS that are potentially suitable for an 
experiment of this nature. Potential candidates for this 
experiment include: 

(1) The Free Electron Laser Driven by the NBS Microtron 

(2) Stanford Superconducting Accelerator/Free Electron 
Laser SCA/FEL 

(3) Lawrence Livermore National Laboratory AT A/ TRW 
Wiggler. 

The NBS CW Microtron Driven FEL (currently under 
development) has the design goals of a normalized transverse 
emittance less than 10 mm-mrad and a momentum spread of less 
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than .03% for operation at 0. 2-3.0 microns. This momentum 
spread is one order of magnitude better than the MKIII IRFEL 
due to the inherent acceleration process of a microtron. 
The hardware and bending magnet configuration downstream of 
the wiggler exit need to be investigated and modeled in 
order to determine a suitable detector location to measure 
the electron beam modulation. 

The Stanford Superconducting Accelerator/Free Electron 
Laser (SCA/FEL) has a normalized transverse emittance of 5 
mm-mrad (10 mm-mrad with a high brightness injector) and a 
momentum spread of approximately .04% for a fundamental 
wavelength of 3u. [Ref. 8] The hardware and bending magnet 
configuration downstream of the wiggler exit also need to be 
investigated and modeled in order to determine suitability 
for a proposed experiment. 

The Advanced Test Accelerator/TRW Wiggler FEL utilizes 
induction linac technology and the laser amplifier FEL 
configuration to amplify a carbon-dioxide laser at 10.6 
microns. This system may be more suitable for detector 
placement due to the absence of optical cavity mirrors and a 
portable beam dump. 

The TRANSPORT equations were generated using a first 
order approximation and an assumed beam spot size of radius 
1 mm at the entrance pole face of the first bending magnet. 
The accuracy of the equations could be improved by utilizing 
the second order fitting capability of TRANSPORT. 
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The system of Figure B-3 was modeled starting at 
location A, assuming a beam spot size of radius 1 mm. A 
more accurate model would incorporate the effects on the 
electron beam envelope within the wiggler. TRANSPORT does 
not have a designated "wiggler" element however, the wiggler 
can be modeled using drift spaces and the "Arbitrary 
Transformation" (type code 14) input option. This model 
should start at the wiggler center where the beam waists x w 
and y w are known. [Ref. 10] 
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APPENDIX A 



BEAM TRANSPORT THEORY 



1. PHASE SPACE 

The transport of a charged particle beam through a 

system of magnets and drift spaces involves the study of a 

fundamental conservation law known as Liouville's theorem: 

. . . under the action of forces which can be derived from 

the Hamiltonian, the motion of a group of particles is 
such that the local density of the representative points 
in the phase space remains everywhere constant. [Ref. 11] 

Liouville's theorem advises the reader that he must look at 
the collective group of particles rather than an individual 
particle in order to access the overall nature of beam 
transport. A single particle moving in a three dimensional 
coordinate system is completely specified if its position 
and momentum components are known, namely x, y, z, P x , Py, 
P z . It is convenient to express this information in a six- 
dimensional space referred to as phase space. A beam 
comprises all individual particles within this phase space, 
neglecting electrostatic repulsion and spin dependent 
effects. The latter two effects would require an 

unmanageable phase space of dimension six times the number 
of particles in the beam plus the additional dimensions 
corresponding to the rotational degrees of freedom of the 
individual particles. Any system of conservative forces 
acting on a beam possesses a Hamiltonian. The macroscopic 
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magnetic forces along the beam trajectory of the proposed 
experiment are conservative thus the beam must obey 
Liouville's theorem. The concept of a six-dimensional phase 
space is a valid and useful utility in the study of beam 
transport since space charge effects become weaker as beam 
energy increase. The six-dimensional phase space 
representing all particles in the beam is considered a 
hypervolume that can change shape but must maintain constant 
volume. The phase space analysis of the beam can be further 
simplified by the choice of the coordinate axis. If the 
three components of motion are mutually independent in real 
space, the motion is confined to the three planes (x, P x ) , 
(y, Py) , and (z, P 2 ) . The six dimensional phase space 
hypervolume has in effect been reduced to three two- 
dimensional phase space areas. A further simplification in 
the analysis is possible by assuming the axial momentum of 
the beam is constant. This is a valid assumption downstream 
of the accelerator sections of the beam trajectory. The 
angular divergence of a particle relative to the beam axis 
is equal to the ratio of the transverse and axial momenta. 
Considering the transverse x-direction the axial divergence 
would be P x /P z or dx/dz. Dimensional phase space (x,P x ) is 
now represented as (x, dx/dz) where dx/dz is the angular 
divergence, a directly observable quantity. Similarly, the 
transverse y-direction phase space representation is 
(Y/dy/dz). The longitudinal phase space is represented as 
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(l,dP z /P z ) where 1 denotes the longitudinal beam extent and 
dP z /P z denotes the fractional momentum spread. [Ref., 11] 

2 . EMITTANCE 

The quality of any electron beam is measured in terms 
of its emittance. Consider a single electron propagating in 
the z-direction as shown in Figure A-l. 



> 


electron having off axis 
< r~ component of velocity 


. i T 




■ ■ - r W c 

electron on beam axis 



Figure A-l Single Electron Motion in the Horizontal Plane 

Emittance is conveniently expressed in terms of the phase 
space parameters x and dx/dz. Using the small angle 
approximation dx/dz can be represented as 0. Electron beams 
produced by accelerators have elliptical phase spaces as 
shown in Figure A-2 where: 

Bjuax ~ the maximum divergence of an individual electron 
x raax = the maximum horizontal extent of a beam electron. 
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Figure A-2 Elliptical Phase Space 

In accordance with Liouville's theorem, phase space is 
conserved. The electrons are within the perimeter of the 
ellipse. Emittance is mathematically defined: 

e = (area of phase space) . (A-l) 

Normalized emittance is defined: 

e n = / (A-2) 

where y is the ratio of the particle energy to the rest mass 
energy and 3 = 1. The area of the elliptical phase space is 
irXmaxBmax therefore: 
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Eft = Xmax^max. 



For example: The Stanford MKIII FEL has the following 

normalized parameters [Ref. 12]: 

horizontal phase space area = 8 nmm-mrad , 
vertical phase space area = 4 tt mm-mrad . 

therefore assuming E = 40MeV: 

e nh = 8 mm-mrad ; e ^ = . 10 mm-mrad , 
z nv = 4 mm-mrad; z v = .05 mm-mrad . 

Assuming elliptical phase spaces, we can assign approximate 
values for Xmax and Ymax and the corresponding values for 
6max and cjjmax. 

Letting Xmax = Ymax = 1mm and energy equal to 40MeV with 
3 = 1 the corresponding values of max and max are 

determined to be: 

Xmax = 1 mm; e max = .10 mrad , 

Ymax = 1 mm; ^ max = .05 mrad . 
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3. 



BEAM WAISTS 



In accordance with electron beam transport design [Ref. 
10], the electron beam is normally focussed at the center of 
the wiggler of the FEL. This yields a beam waist in the 
center of the wiggler as shown in Figure A-3. 




Figure A-3 Beam Focus at Center of Wiggler 



where 



x w = horizontal beam waist 

y w = vertical beam waist 

Iv = length of wiggler (1 meter) 

The horizontal beam waist size is related to the length of 
the wiggler and the horizontal emittance by the relation: 

x w = kw e h/2 (A.-3) 
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The vertical beam waist size is related to the vertical 



emittance and the intrinsic betafunction, 3 w , by 
relation: 



Yw = ^v^w 



(A 



where 



the 



4) 



6 w = electron radius of curvature at 
center of wiggler. 



APPENDIX B 



TRANSPORT-A COMPUTER PROGRAM FOR DESIGNING 
CHARGED PARTICLE BEAM SYSTEMS 

1. INTRODUCTION TO TRANSPORT 

The study of static-magnetic beam transport systems led 
to the development of the original first-order TRANSPORT 
computer program. TRANSPORT was written in BALGOL at 
Stanford Linear Accelerator Center (SLAC) in 1963. It has 
been refined in the succeeding years by contributions from 
various laboratories throughout the world. Subsequent 
improvements include the translation into FORTRAN and a 
second-order fitting capability. The version of TRANSPORT 
used in this thesis is the PC version dated May, 1985 from 
D. Carey, Fermi National Accelerator Laboratory (FNAL) , 
compiled with the IBM version of professional fortran by 
Ryan-McFarland. [Ref. 13] 

TRANSPORT utilizes matrix theory methods to step 
sequentially through a beam line system consisting of 
magnets, drift spaces, and special magnet configurations. 
The beam line system and the input beam parameters are 
specified by the user and an output beam matrix and overall 
system transfer matrix are generated. Each individual step 
through the system is concidered an element. Each element 
is given by a sequence of items separated by spaces and 
terminated by a semicolon. These items are a type code 
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number, a vary field, the physical parameter and optional 
labels. The type code number identifies the element (i.e., 
magnet, drift space, pole face rotation) . A vary field 
allows the user to adjust the physical parameters of an 
element if the problem involves fitting. Physical 
parameters are quantities which describe the physical 
element. Optional labels assist in identifying beam 
elements and parameters in the output. A summary of 
TRANSPORT type codes is shown on page 77. [Ref. 13] 

The initial objective in writing the input program will 
be to specify the beam line, establishing the central 
trajectory. The system to model by TRANSPORT is shown in 
Figure B-l. This system consists of two bending magnets and 
five drift lengths where 

L 0 = initial drift length , 

L = effective pathlength through bending magnet , 
L]_ = drift length 1 , 

L 2 = drift length 2 , 
a = bending angle . 

2. SAMPLE PROGRAM DEVELOPMENT [Ref. 13] 

Figure B-2 is a sample input program that establishes 
the beam line of Figure B-l, assuming the initial drift 
space, L 0 equals zero. Refer to Figure B-2 for the 
remainder of this section. 
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Figure B-l System Configuration 



(1) ’STANFORD MKIII IRFEL : REVISED 10/16/87’ 

0 



(2) (CARTESIAN COORDINATE SYSTEM; x.y.s) 



13. 


’COOR' 


12. ; 






(SET 


UN JTS/scaled 


units ) 






15. 


’ xymm’ 


1 . 


’ mm ' 


0. 1 


15. 


’ EM^V ’ 


11 . 


•MeV 


0.001 


15. 


’ long ' 


8. 


’ mm ' 


0.001 


15. 


' leng * 


5. 


’ mm ’ 


0. 1 



(4) (BEAM PARAMETERS AT FIRST BENDING MAGNET POLE-FACE) 

( NORMAL I ZED BEAM EMITTANCE: Enh=8min-mrad , Env=4mm-mrad ) 
1.0 ’BEAM’ 1.0 0.1 1.0 0.05 .0005 0.0 40.0 

(monoenergetic beam) 

(5) 16. 5. 5.0 ; (half-gap). 

(default value of Kl=.5) 

(6) (FIRST BENDING MAGNET) 

2 . 0 . ; 

4. ’ FNDJ ’ 38.347 6.33175 0. ; 

2. 11.25 ; 

13. 1. ; (PRINT BEAM MATRIX) 

13. 4. ; (PRINT XFORM MATRIX) 

(7) (DRIFT SPACE rtl) 

3. ’ DFT1 ’ 327.03 ; 

13. 1 . ; 

13. 4 . ; 



(std coord rot w/11.25 bend) 
6.83175 0. ; 



(9) (DRIFT SPACE 112) 

3. ’ PFT2 ' 95.30 

13 . 1 . ; 

13 . 4 . ; 

(10) SENTINEL ; 



(8) (SECOND BENDING MAGNET) 



20. 


180. 


’LEFT’ 


2. 


11.25 


t 


4 . 


’ BND2 ’ 


38. 347 


n 

L. . 


0. ; 




2n . 


-130. 


; 


13. 


1 . ; 




13. 


4 . ; 





Figure B-2 Sample Input Program 
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The individual elements and there associated type codes will 
be explained in the succeeding paragraphs. 

(1) 'STANFORD MKIII PEL: revised 10/16/87' is the 
arbitrary title of the sample input program. The zero 
immediately below the title indicates the start of a new 
problem. 

(2) A cartesian coordinate system was chosen for the 
proposed experiment design problem. The coordinate system 
is established using the type code 13. The statement [13. 
12. ;] will yield a beam line reference trajectory that will 
start at the origin and proceed along the positive z-axis. 
The y-axis will point up and the x-axis will point to the 
left. Statements enclosed by ( ) are comment statements to 
assist in following the code. 

(3) The scaling of units is accomplished via type code 
15. All lengths in Figure B-l are to be specified in 
millimeters and energy is to be specified in MeV. TRANSPORT 
has the standard set of units as shown on page 78. 
Conversion factors are applied to code digits 1, 11, 8, and 
5 to establish an appropriate unit system. A summary of 
TRANSPORT type codes is found on page 77 . 

(4) The input beam, denoted by the type code 1, is 
specified in terms of its phase space and average momentum 
parameters. The phase space input is given in terms of the 
semi-axis of a six-dimensional beam ellipsoid representing 
the phase space variables x, 0 , y, $ , 1, and 6 where 
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X 



= one-half the horizontal beam extent (1.0mm) , 

6 = one-half the horizontal beam divergence (O.lmrad) , 

y = one-half the vertical beam extent (1.0mm) , 

<J> = one-half the verticle beam divergence (.05mrad) , 

1 = one half the longitudinal beam extent (.0005mm) , 

6 = one-half the momentum spread (0.0%) , 

P = average momentum of the beam (40MeV/c) . 

NOTE: The above units were scaled using type code 15. 

(5) and (6) Bending magnet parameters are specified in 
terms of the wedge magnet, type code 4; pole-face rotation, 
type code 2; and a gap setting, type code 16. A wedge 
bending magnet implies that the central trajectory of the 
beam enters and exits perpendicularly to the pole face 
boundaries. A wedge magnet can be modified into a square 
edge magnet using the pole face rotation element on either 
the entrance or exit face of the magnet. The field 
boundaries for bending magnets are shown on page 73. To 
modify a wedge bending magnet into a square for a bend to 



the right. 


Pi 


is set to 


zero and 


B 2 is set 


to 


11.25 


degrees . 


The 


pole face 


rotation 


type code 


2 . 


must 



immediately precede or immediately follow the bending magnet 
(type code 4.0). There are four first order parameters to 
be specified for the wedge magnet; the type code, the 
effective length of the central trajectory, the central 
field strength and the field gradient. The bending magnet 
equations and parameters for the system are listed below: 
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L = p^a where 

L = effective pathlength through magnet , 
p 0 = bending radius; (195.30 mm) , 
a = bend angle; (11.25 degrees) , 

L = (195. 30mm) (11. 25deg) ( rad/180deg) = 38.347 mm . 

(NOTE; To simulate a gradual variation of the fields at 
the ends of the bending magnets a fringing field 
element is provided.) 

B ( 0 ) = 3 3 . 3 5 6 ( P/ p 0 ) where 

B ( 0 ) = central field strength , 

P = central momentum; (40 MeV/c) , 

B ( 0 ) = 33.356 (40MeV/c) /195 . 30mm = 6.83175 kG , 

n = magnetic field gradient (n = 0 assumed) . 

Code [13. 1. ;] yields the printing of the current beam 
( o) matrix. Code [13. 4. ;] yields the printing of the 
current transformation matrix, Rl. These codes will be 
addressed in the matrix formulation of TRANSPORT. 

(7) and (9) A drift length (type code 3.0) is a field 
free region through which the beam passes. It is denoted by 
the type code and the effective drift length; 

Li = drift length 1 (327.03 mm) , 

L 2 = drift length 2 (95.30 mm) 

(8) The second bending magnet differs from the first 

in that is set to 11.25 degrees and i s set to zero 

degrees. In addition, a coordinate rotation element (type 
code 20.) is required to rotate the transverse coordinates x 
and y through an angle about the z-axis. A bend to the left 
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(looking in the direction of beam travel) is accomplished by 
rotating the x,y coordinates by 180 degrees. 

(10) The word "sentinel" denotes the end of the input 
program. 



3. MATRIX FORMULATION OF TRANSPORT [Ref. 13] 



A. Phase Space Vector 

At any specified position in the system an arbitrary 
charged particle is represented by a vector X, whose 
components are the positions, angles, and momentum of the 
particle with repect to the central trajectory. The X 
vector is represented: 

’ x ‘ 



X = 



0 

Y 

4 > 

1 

6 



where 

x = the horizontal displacement of the arbitrary ray 
with respect to the assumed central traj ectory , 

0 = the angle this ray makes in the horizontal plane 

with respect to the assumed central trajectory , 

y = the vertical displacement of the ray with respect to 
the assumed central trajectory , 

<p = the vertical angle of the ray with respect to the 
central trajectory , 

1 = the path length difference between the arbitrary ray 

and the central trajectory , 

6 = Ap/P; the fractional momentum deviation of the ray 
from the assumed central trajectory. 
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B. The Transfer Matrix R 

The passage of a charged particle through an element 
(magnet, drift space, rotation, etc.) may be represented by 
Equation B-l: 



X(l) = RX ( 0 ) 



(B-l) 



where 

X ( 0 ) is the initial coordinate vector, 

X(l) is the final coordinate vector, 

R is the transformation matrix of the element. 

In a system of several magnets, pole face rotations, and 
drift spaces, R is replaced by the product matrix R(t) : 

R (t) = R(n) . . .R(3)R(2)R(1) . (B-2) 

R(t) is the product of the individual element matrices of 
the system. Equation B-3 represents the matrix equation for 
a system: 



Xi(t) 



6 

I RijXj (0) 



(B-3) 



where 



Xi = X, 

x 2 = e. 
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x 3 = y, 



x 4 = <P, 

X5 = 1 , and 
X 6 = 5 . 

This product matrix is automatically calculated by the 
program and is called TRANSF 0 RM 1 . TRANS FORM 1 can be printed 
when desired using the statement [ 13 . 4 . ;] (e.g., paragraph 

(6)). For static magnetic systems possessing midplane 
symmetry, the six simultaneous linear equations represented 
by Equation B -3 may be expanded in matrix form as [Ref. 14 ]: 



x(t) 




R 11 


R 12 


0 


0 


0 


R 16 




x o 


e(t) 




R 2 1 


R 22 


0 


0 


0 


R 2 6 




e 0 


y(t) 




0 


0 


R 3 3 


R 34 


0 


0 




y 0 


4 >(t) 




0 


0 


R 4 3 


R 44 


0 


0 




^0 


i(t) 




R 51 


r 52 


0 


0 


1 


R 56 




lo 


6(t) 

- 




0 


0 


0 


0 


0 


1 




- 6 °_ 



The transformation is from an initial position t = 0 
to a final position x = t measured along the assumed central 
trajectory. The final coordinate vector X(t) is the matrix 
product of the individual coordinate vector and the overall 
system transformation matrix. The individual Rj_j elements 
of the matrix are dimensionally consistent with the input 
units specified in the transport code. The individual 
transformation matrices needed to represent the system of 
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Figure B-l are shown on pages 74-76. It is desired to 

determine and plot x(t) vs l(t) for the overall system. 
Since the system is symmetric, x(t) and l(t) have no 
dependence upon y Q or 4> 0 , thus the general form of the 
tranf ormation matrix is reduced to: 



x(t) 




R 11 r 12 0 r 16 




x o 


0(t) 




r 21 r 22 0 r 2 6 




e o 


Kt) 




r 51 r 52 1 r 56 




lo 


6(t) 




0 0 0 1 




V 



The wedge bending magnet is reduced to: 



cos a Po sin a 0 

-hsin a cos a 0 

-sin a -p 0 (l-cosa) 1 

0 0 0 



P Q (1-cosa ) 
sin a 

- p Q ( a-sin a) 
1 



The pole face rotation matrix is reduced to: 



1 0 

(tana)/ Po 1 

0 0 

0 0 



0 

0 

1 

0 



0 

0 

0 

1 
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The drift space matrix is reduced to: 



1 

0 

0 

0 



L 

1 

0 

0 



0 

0 

1 

0 



0 

0 

0 



It is desired to determine the X matrix at various locations 
along the central trajectory. Figure B-3 denotes the system 
locations to be evaluated: 




Figure B-3 System Configuration Evaluation Points 
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The matrix representation of X B (t) is the product of the 
pole face rotation matrix, wedge bending magnet matrix, and 
the initial coordinate vector: 



x B (t) 




i 


0 


0 


0 




cos a p 0 sina 0 


P 0 ( 1-cosa) | 


x ° 


0 B (t) 




tan cx/ pQ 


1 


0 


0 




-hsina cos a 0 


sina 


e o 


l B (t) 




0 


0 


1 


0 




-sin a “p 0 ( 1-cos a) 1 


- p Q ( a-sina)| 


loi 


1 

Cb 

Gd 

rt 

i . 




0 


0 


0 


1 

- 




0 0 0' 


1 

J 


6 o- 



(pole face rotation) (wedge bending magnet) 



The matrix product reduces to: 



x B (t) 




cos a 


p Q sina 0 


p Q ( 1-cos a) 




— — t 

x o 


0 B (t) 




0 


1/cos a 0 


tan a 




e o 


l B (t) 




-sin a 


“ Po ( 1-cos a) 1 


- ( a-sin a) 




-*-0 


6 B (t) 




0 


0 0 


1 




6 

o 



The matrix representation of X c (t) is X B (t) multiplied by 
the first drift space matrix: 

















- 




n 


- 


x c (t) 




1 


L 1 


0 


0 




cos a 


p^sin a 0 


( l-cos a) 




x o 


0 cCt) 




0 


1 


0 


0 




0 


1/cos a 0 


tan a 




0 O 


l C (t) 




0 


0 


1 


0 




-sim 


-(^(1-cosa) 1 


- p 0 ( a-sin a) 




lo 


<5 C (t) 




0 


0 


0 


1 




0 


0 0 


1 




6 

o 
















- 








_J 



(first drift) X B (t) 
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The matrix product reduces to: 



x c (t) 




cos a 


p 0 sina+L 1 /cosa 


0 


p Q ( 1-cosa) +L^tan a 




x o 


CD 

o 

r+ 




0 


1/cosa 


0 


tan a 






6 o 


i c (t) 




-sin a 


- p Q ( 1-cosa) 


1 


- p D ( a-sina) 






lo 


<5 c (t) 

J 


1 


0 


0 


0 


1 






_ 6 oj 


The matrix representation of 


X D (t) 


includes a 


coordinate 



rotation which will in effect bend the beam to the left. In 
matrix form this is represented by changing the bending 
radius and bending angle to its negative value (i.e., p Q — > 
-p D / a -»“a) • The wedge bending matrix for a bend to the left 
is modified as shown : 



cos a P Q sin a 0 

-hsina cos a 0 

sin a po(l-cosa) 1 

0 0 0 



p Q ( 1-cosa) 
-sina 

-p 0 ( a-sina) 
1 



The pole face rotation matrix retains its identical form 
because tan (-a) = -tan (a)- Multiplying Xc(t) by the pole 

face rotation matrix and the wedge bending magnet matrix 
yields : 
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X D (t) 




1 


2 Po tan cd'Li/ cos 2 a 


0 


Zl 


1 1 


0D( t ) 




0 


1 


0 


0 


1 


l D (t) 




0 


Z2 


1 


Z3 




<S D (t) 




0 


0 


0 


1 














J 





where : 

Z 1 = -1 p Q + 2p 0 /cosa + L]_sina/cos 2 a , 
Z2 = Zl , 



Z3 = -2p 0 a + 2 p Q tana + L]_tan 2 a . 



The matrix representation of X]?(t) is X[)(t) multiplied by 
the second drift space matrix: 



x E (t) 

0 E (t) 




“ 

1 Z4 0 

0 1 0 


O tN) 

U1 

I 




l 

CD X 

0 0 

l 


!e ( t ) 

' E 




0 Z 6 1 

0 0 0 

L 


Zl 

l 




lo 

1 ^ o 



where: 

Z4 = 2 p Q tan a + L 1 /cos 2 a + l 2 , 

Z5 = -2 p Q + 2 p 0 /cos a + L^sin a/cos 2 a , 
Z 6 = Z 5 , 

Z7 = -2 p Q a + 2 P Q tan a + L]_tan 2 a . 
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In equation form, the matrix products are represented as: 



x B (t) = (cos a) x Q + ( Po sin a) ^ + p Q ( 1-cos a) 5 0 (M-l) 

OgCt) = (1/cos a) 0 O + (tana) 6q (M-2) 

l B (t) = (-sin a) x Q - ^(l-cosaJGo + 1 0 - Pb( a - sina)<5 0 (M-3) 
6 B (t) = 6 0 (M-4 ) 



x c (t) = (cos a )x 0 + ( Po sin a + Lx/cos a) Qq + [^(l-cosa) (M-5) 

+ L x tan a] <5 q 

^(t) = (1/cos a) 0 O + (tana) 6 q (M-6) 

lc ( t) = (-sina)x Q - p Q ( 1-cos a) 0 O + lo - P Q (a-sina) 6 Q (M-7) 

6 c (t) = 5 0 (M-8) 

x c (t) = x B (t) + I^G^t) (M-5A) 

0 c(t) = %(t) (M-6A) 

l c (t) = l B (t) (M-7A) 

6 cCt) = <5 B( fc ) (M-8A) 

^D(t) = *o + ( 2 Qptan a + L : /cos 2 a) ^ + (-2^ + 2 p 0 /cos* (M-9 ) 

+ L]_sin c/cos 2 a) <5^ 

e D( t ) = 6 0 (M-10) 

lD( t ) = (“2 to + 2 to /cos a + I^sin c/cos 2 a)9 0 + 1 Q (M-ll) 

+ (-2 pQa + 2p 0 tana + L 1 tan 2 a)6 Q 

6 D( fc ) = 6 o (M-12) 
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X D (t) = 


( l/cosa) x c ( t) + (p Q sin ) ^(t) - p Q (l- 


cos a) 6c (t) 

(M-9 A) 


II 

Q 

CD 


(cos a) 0 C (t) - (sina) 6c( t ) 


(M-10A) 


l D (t) = 


(tan a) x c (t) + ( l-cosa)0 c ( fc ) + 1 C 

- p 0 (a “ sin ) 6c(t) 


(M-11A) 


II 

Q 

to 


O 

ft 


(M-12A) 


x E (t) - 


X Q + (2p 0 tana + L-,/cos 2 a + L 2 ) e 0 

+ (-2p 0 + 2p 0 /cos a + L^sina/cos^ot 


(M-13 ) 

) 6 o 


0 E (t) = 


9 o 


(M-14 ) 


l E (t) = 


( -2 p 0 + 2 q-)/ cos a + Lxsina/cos 2 a) 0 p + 1 Q 
+ (-2p 0 a + 2p 0 tana +L]_tan 2 a) o Q 


(M-15 ) 


) = 


6 o 


(M-16) 


x E (t) = 


x D (t) + L 2 e D (t) 


(M-13A) 


II 

W 

CD 


^(t) 


(M-14A) 


II 

4^ 

W 

rH 


in(t) 


(M-15A) 


6 E (t) = 


6 D( fc ) 


(M-16A) 


c. 


Beam Matrix 





The phase ellipse beam (sigma) matrix is used to 
inform the user of the beam parameters at any location in 
the system. The beam matrix will represent the bundle of 
particles or rays constituting the phase space of the input 
beam at desired locations in the system. Particles of the 
beam are assumed to lie within the boundaries of the 
ellipsoid with each point within the ellipsoid representing 
a possible ray. The equation of an n-dimensional ellipsoid 
may be written in matrix form as: 



65 



(B-4) 



X ( 0 ) T a (0)" 1 X(0) = 1 , 



where: 



X(0) T = transpose of coordinate vector X(0), 

a(0) = a real positive definite, symmetric matrix. 

As a particle passes through the system, it 
undergoes the matrix transformation of Equation B-l. 
Combining this transformation with the equation of the 
initial ellipsoid and using RR -1 = I : 

X(0) t (R t R t-1 ) a(0) -1 (R -1 R)X(0) = 1 . (B-5) 

Equation B-5 can be rewritten as: 

[RX(O) ] T [Ra(0)R T ]“ 1 [RX(0) ] = 1 . (B-6) 

The equation of the new ellipsoid after the transformation 
becomes : 



X(l) T a(i)-lx(l) = 1 , 



(B-7) 



where: 



a ( 1 ) = Ra ( o ) R t 



(B-8) 
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Thus the beam may be represented at any point in the system 
by Equation B-8. Note: a (0) is the input beam phase space 
specified by the user. The beam matrix is printed using the 
statement [13. 1. ;]. The projection of the semi-axis of 
the ellipsoid upon each of its six coordinate axes is 
printed in a vertical array, and the correlations among 
these components are printed in a triangular array. The 
phase ellipse beam matrix in the printed output form is 
shown in Figure B-4. 









X 


0 


y 




l 


X 


✓o(ll) 


CM 












0 


/a(22) 


MR 


r(21) 










y 


^(33) 


CM 


r(3l) 


r(32) 










✓0(44) 


MR 


r(41) 


r (42) 


r (43) 






i 


^0(55) 


CM 


r (51) 


r (52) 


r (53) 


r(54) 




6 


/c( 66) 


PC 


r(61) 


r (62) 


r(63) 


r(64) 


r(65) 



where: r(ij) » -Hiiil — . 

[c(ii)a(j j)] * 

As a result of the fact that the 0 matrix is positive definite, the 
r(ij) satisfy the relation 

|r(ij)| £ 1 . 



Figure B-4 Beam matrix 



67 



(Note: units will be labeled as specified by type code 15.) 

The individual sigma elements are defined: 

Jo (11) = x max = the maximum half-width of the beam 

envelope in the x (bend) -plane at the 
point of the printout, 

Jo (22) = ©max = the maximum half-angular divergence of 

the beam envelope in the x(bend) plane, 

Jo (33) = y max = the maximum half -height of the beam 

envelope, 

Jo (44) = (fonax = the maximum half-angular divergence of 

the beam envelope in the y (non-bend) 
plane, 

Jo (55) = l max = one-half the longitudinal extent of the 

bunch of particles, 

Jo (66) = 6 = the half-width 1/2(AP/P) of the momentum 

interval being transmitted by the system. 



Note: the units appearing next to the Jo (ii) in the 

TRANSPORT printout are the units chosen for coordinates x, 
0 Yi 4> / 1/ and 6 respectively. 



An example of the (x,6) plane ellipse is illustrated 
in Figure B-5. 

Consider a two dimensional (x, 0 ) plane projection 
of the general 6-dimensional ellipsoid. Let: 



an a 21 



a _ 



°21 



°22 
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Figure B-5 A Two Dimensional Beam Phase Ellipse 



be a real, positive definite, symmetric matrix. The inverse 
of the o -matrix exists and is defined by: 



a 



1 



1/e 2 



022 ~°21 

” °21 a ll 



where e 2 is the determinant of 

The two-dimensional coordinate vector and its 
transpose are represented by: 

x 



X 



e 



x T 



(x, 0) . 
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The expansion of the matrix equation 



X T a” 1 X = 1 

is the equation of an ellipse, therefore 

°22 x2 “ 2 a 21 X 0 + o^e 2 = e2 = deta . 

The area of the ellipse is given by: 

1/2 

A = tt (deta ) = Trx max e j_ n t = TrXi n t®max = i;e . 

The correlation between x and 0 (the orientation of the 
ellipse) depends on the off diagonal term a 21* This 
correlation is defined as: 

a 21 

r 21 = r 12 = — 

' /a ll a 22 

The correlation, r, measures the tilt of the ellipse and the 
intersection of the ellipse with the coordinate axis. The 
beam ellipse in the (x, 0) plane will be valuable in the 
analysis of the beam divergence at various points in the 
system and also in the fact that it reveals the change in 
the longitudinal extent of the beam via element /a (55) . 
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4. INSTALLATION OF TRANSPORT ON THE IBM PC [Ref. 15] 

The following requirements and procedures are necessary 
for the loading and running of TRANSPORT: 

Requirements : 

- config.sys file must have files >= 15; buffers >= 15. 

- disk set A containing .EXE and .BAT files in DOS/BACKUP 

- 64 OK memory 

- hard disk; lOM-byte or larger 
Procedures : 

- move to the C directory: C: [ENTER] 

- create a TEMPO subdirectory: MD\TEMPO [ENTER] 

- move to the TEMPO subdirectory: CD\TEMPO[ ENTER] 

- RESTORE .EXE and .BAT files: RESTORE A: C: [ENTER] 

(insert disk #1 and #2 as directed) 

- move to the C directory: CD\C: [ENTER] 

- create a TRANS subdirectory: MD\TRANS [ENTER] 

- move to the TEMPO subdirectory: CD\TEMPO [ENTER] 

- copy files from TEMPO to TRANS: copy *.* \TRANS [ENTER] 

- delete files from TEMPO subdirectory: del *.* [ENTER] 

- move to the C directory: CD\C: [ENTER] 

- remove TEMPO subdirectory: RD\TEMPO [ENTER] 

Limitations : 

- the PC version of TRANSPORT will allow for 1000 elements 
and 4500 data values 

- accuracy is four significant figures 
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To run TRANSPORT: 



*- TRANSGO.BAT is the command file to run TRANS.EXE 

- type the following: 

TRANS GO FILE1 FILE2 FILE3 
FILE1: input data file name (users file name) 

FILE2: output print file name (users file name) 
FILE3: output summary file name (users file name) 

- example: TRANSGO MK3FELA MK3FELA. OUT MK3FELA.JNK[ ENTER] 

* Note: TRANSPORT is currently loaded in the IBM PCs in the 
Physics Library and Physics Department Office under 
the subdirectory C:\TRANS. X. Maruyama currently 
holds the TRANSPORT software on 5-1/4 disks. 
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First-order wedge bending magnet matrix 
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path length of the central trajectory. 



TRANSPORT calculation is as follows 




75 




76 



15. 8. ’ ’ ; type code entry (if used) preceding the BEAM (type code 1 
card. If no 15. 8. entry is made, the units of L. will automatically be 
in metres (standard TRANSPORT units) . 



Table l: Summary o( TRANSPORT type codes 
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The various units that may be changed are: 



Code 

Digit 


Quantity 


Standard 

TRANSPORT Unit 


Symbols used in 
S LAC-75 


1.0 


horizontal and vertical trans- 
verse dimensions, magnet 
apertures, and misalignment 
displacements 


cm 


x,y 


2.0 


horizontal and vertical angles 
and misalignment rotation 
angles 


mr 


e,<p 


3. 0 


* \ 

vertical beam extent (only) ’ 
and bending magnet gap height 


cm 


y 


4.0 


* ) 

vertical beam divergence ' 
(only) 


mr 


<P 


5. 0 


pulsed beam length and wave 
length in accelerator 


cm 


i 


6. 0 


momentum spread 


percent (PC) 


6 


7.0 


bend, pole face rotation, and 
coordinate layout angles 


degrees (DEG) 




8.0 


length (longitudinal) of ele- 
ments, layout coordinates, 
and bending magnet pole face 
cur vatu res 


metres (M) 


t 


9. 0 


magnetic fields 


kG 


B 


10. 0 


mass 


electron mass 


m 


11. 0 


momentum and 


GeV/c 


P(0) 




energy gain in accelerator 
section 


GeV 


AE 



*) These codes should pot be used if the coordinate rotation (20.0) type code is 
used anywhere in the system. 
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